Resonant Two-Magnon Raman Scattering in Cuprate Antiferromagnetic Insulators 
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We present results of low-temperature two-magnon resonance Raman excitation profile measure- 
ments for single layer Sr2Cu02Cl2 and bilayer YBa2Cu306+i antiferromagnets over the excitation 
region from 1.65 to 3.05 eV. These data reveal composite structure of the two-magnon line shape 
and strong nonmonotic dependence of the scattering intensity on excitation energy. We analyze 
these data using the triple resonance theory of Chubukov and Frenkel (Phys. Rev. Lett., 74, 3057 
(1995)) and deduce information about magnetic interaction and band parameters in these materials. 



PACS numbers: 78.30.Hv, 75.50.Ee, 74.72.Bk 

The insulating phases of the high- Tc materials are an- 
tiferromagnets (AF) with spin S = i, localized on Cu 
atoms of the Cu02-planes, characterized by superex- 
change constant J « 120 meV and a Neel temperature 
Tn « 300 K The magnetic interactions in the 

Cu02-planes can be described by a Hcisenberg Hamilto- 
nian on the 2D square lattice H — JJ2<i j>i^i '^j ~ j)^ 
where Si is the spin on site i and the summation is over 
nearest-neighbor Cu pairs. 

Two-magnon (2M) Raman scattering (RS) probes 
mainly short wavelength magnetic excitations in the AF 
lattice. The scattering process involves a virtual charge- 
transfer excitation. The traditional Hamiltonian for de- 
scribing the interaction of light with spin degrees of 
freedom is the Loudon-Fleury Hamiltonian = 
'^J2<i,j>i^i ■ • • Sj, where and are 

the polarization vectors of incoming and outgoing pho- 
tons, a is the coupling constant and Rij is a vector con- 
necting two nearest-neighbor sites i and j. Shastry and 
Shraiman 0| have derived this Hamiltonian in the large- 
ly Hubbard model approach for non-resonant excitation 
energies. Including final state interactions in this 2M 
Raman scattering theory results in a peak with Raman 
frequency shift w near 2.7 — 2.8J in Big scattering ge- 
ometry for a single layer 2D lattice of D^h. symmetry 
|^|-^. This allows one to determine the exchange in- 
teraction constant J directly from the peak energy ob- 
served in the Raman spectrum. In fact, the Raman data 
from the Big 2M scattering peak yielded the first es- 
timate of the exchange interaction constant in La2Cu04 
|-|ll|, Gd2Cu04 H, RBa2Cu306+5 (R-Y Eu @, 
Pr and Nd |5|), Sr2Cu02Cl2, LaGdCu04, Nd2Cu04, 
(Ca,Sr)Cu02 |11|]. The 2M peak has also been observed 
in the electron-doped materials ^6 17 1 and in under- 
doped YBa2Cu306+5 and YBa2Cu408 superconductors 
The temperature dependence of the 2M scattering 
has been studied in Ref. |l4[ . 

The dependence of 2M scattering intensity on the in- 
coming excitation energy tOi has been discussed in a num- 
ber of earlier works ||lq-|2q| . Comparison of optical con- 



ductivity data with these Raman measurements shows 
that the maximum strength of the 2M scattering occurs 
for the excitation energies substantially higher than the 
optical band edge that marks the onset of particle-hole 
excitation across the charge transfer gap (CTG). 

In this Letter we report low-temperature two-magnon 
resonance Raman excitation profile (RREP) measure- 
ments for single Cu02 layer Sr2Cu02Cl2 and bilayer 
YBa2Cu306-i-5 AF insulators, covering the broad exci- 
tation region from 1.65 to 3.05 eV. These data reveal 
composite structure of the 2M line shape and strong non- 
monotonic dependence of the shape and 2M scattering 
intensity on the excitation energy uji. We analyze these 
results within the framework of the Chubukov-Frenkel 
approach |8[|2^, show that triple resonance is the lead- 
ing process for resonant magnetic Raman scattering, and 
extract information about the charge transfer band and 
magnetic interactions in these cuprates. 

The investigations were performed on single crystals 
grown as described in Refs. p^ , p9t . Spectra reported here 
were excited with different Ar^ , Kr'^ and dye laser lines 
in x'y', x' x' , xy and xx pseudo-back scattering geome- 
tries from an [001] crystal face. With assumed D^h sym- 
metry, this gives Raman spectra of Big + A2g, Aig + B2g, 
B2g + A2g and Big H- Aig symmetry respectively. The 
A2g component in the 2M excitation region is known to 
be quite weak p^ ]. The single crystals were mounted 
in a continuous helium flow optical cryostat (Oxford In- 
struments) and kept at 5 K. Less than 15 mW of the 
incident laser power was focused in a 50 /im diameter 
spot on the a6-plane crystal surface. Light collected in 
pseudobackscattering geometry was focused with achro- 
matic optics (1:3 magnification) into the 0.5 mm (un- 
derfilled) entrance slit of a triple grating spectrometer 
(SPEX 1877). 150 lines/mm gratings were used to cover 
a broad spectral range. The dispersed light was detected 
by a liquid nitrogen cooled 512 x 512 pixel CCD cam- 
era (Photometries). Spectral corrections to the detector 
sensitivity and the frequency dependence of the collec- 
tion optics and the spectrometer system were made using 
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a calibrated uniform light source sphere (USS-600 from 
Labsphere). We present the spectra in relative units of 
Raman cross section per unit solid angle per unit volume 
per unit frequency shift. The scattered intensity (//) 
is corrected to obtain the cross section via cPa/dujdVl (x 
If{ai+af)n^/{TiTf{Ii/uJi)). This expression corrects for 
transmission coefficients at the sample-helium interface 
for the incident (T^) and scattered (T/) photons, absorp- 
tion of the incident (a^) and reabsorption of the scat- 
tered (a/) photons and also the scattered photon solid 
angle correction at the sample-helium interface. These 
optical parameters were acquired from polarized reflec- 
tivity and ellipsometry measurements. We have used the 
80 K optical data from Refs. for YBaaCugOe+s 

and Sr2Cu02Cl2 crystals respectively. The maximum ef- 
fect of these corrections in the energy range 1.5 to 3 eV 
was ~ 300%. 

The Raman spectra excited with photon energies be- 
tween 1.9 and 3.1 eV in x'x' and x'y' geometries and 
dependence of the Big 2M peak intensity on the excita- 
tion energy are shown in Figs. |l| - 

The features of primary interest are the following: 

(a) Resonant magnetic RS gives the main contribution 
in the Big channel. However, as first shown by Sulewski 
et al. |1^, it contributes to the Aig and B2g scatter- 
ing geometries as well. For these geometries RS bands 
are broad and have less structure. The Aig intensity is 
weaker than in the Big channel, and the B2g intensity is 
weaker than Aig. 

(b) In the Big channel the 2M band is composite and 
contains two asymmetric lineshapes. Sr2Cu02Cl2 shows 
a band with peaks at 2970 and 3900 cm"-'^, where the sec- 
ond peak is well resolved for certain excitation energies. 
For YBa2Cu306+5 the 2M band is peaked at 2720 cm-\ 
and the second peak at 3660 cm~^ is resolved only for the 
lower excitation energies. The 2M band becomes much 
broader and more symmetric when the excitation energy 
approaches 3 eV. 

(c) In the Aig channel for Sr2Cu02Cl2 the 2M band 
intensity at lower frequencies is a combination of a linear 
and a cubic power law. For YBa2Cu306+5 the 2M band 
intensity has flat low-frequency tail and a threshold at 
about 1750 cm"-'^. The sharp peaks at the low energy 
tails are due to resonant multi-phonon scattering that 
becomes strongly enhanced for excitations closed to the 
CTG energy. 

(d) The 2M scattering intensity strongly depends on 
the incoming photon energy. For YBa2Cu306+5 the Big 
peak intensity at 2720 cm^^ is enhanced by three orders 
of magnitude when the excitation approaches 3 eV, that 
is about 1.3 eV higher than the first peak in £2 (See Fig. 

, and a weaker resonance is observed at about 2.1 eV, ~ 
0.32 eV higher than the peak in €2- For Sr2Cu02Cl2 the 
RREP at 2970 cm~^ shows resonance at 2.4 eV, about 
0.44 eV higher than the €2 peak at 1.95 eV, and a steep 
rise to the second strong resonance in the UV region, 



outside the accessibility of our experimental technique. 

These results can be understood in terms of the 2M 
resonant RS theory proposed by Chubukov and Frenkel 
[p|,p7t that suggests a leading resonant RS diagram that 
is different from those leading to the Loudon-Fleury the- 
ory. The proposed process is the following: the incoming 
photon LOi creates a virtual electron-hole state consist- 
ing of an excitation across the CTG, then the fermions 
emit or absorb two magnons (one by each fermion) with 
momenta q and -q before recombining by emitting an 
outgoing photon with the energy w/ = coi — co. The inter- 
nal frequency integral of the corresponding diagram has 
three denominators (See Fig. 2(b) and Eq. (5) in Ref. 
p7| ) that could vanish simultaneously for a certain region 
of Ui and ujf. If the conditions for simultaneous vanishing 
of all three denominators are satisfied, this process gives a 
resonant RS enhancement known as triple resonance . 
The triple resonance enhancement conditions depend on 
the dispersion of the fermionic bands and of the magnons. 
They have been calculated in Ref. for the one-band 
Hubbard model within the spin density wave formalism. 
That model description of the electronic states at half 
filling predicts no triple resonance conditions for excita- 
tion frequencies uji around the CTG energy 2A. The 
enhancement for Raman shifts cu between ^ 2A and 4 J 
occurs twice: first when uji lies between ~ 2.4 and 5.6J 
above 2A and second, even stronger, when uJi approaches 
the upper edge of the electron-hole excitation band (See 
Fig. 3 in Ref. ||). 

This triple resonance process can explain each of the 
observed features (a)-(d): 

(a) The triple resonance contributes to all three exper- 
imentally observed scattering symmetries: Big, Aig and 

B2g- 

(b) The largest contribution from triple resonance scat- 
tering is in the Big channel, where, due to the geomet- 
ric form factor for magnon emission, mainly short wave- 
length magnons from the vicinity of the magnetic Bril- 
louin zone boundary participate in the process. The 2M 
density of states from the vicinity of these points diverges 
at the maximum 2M energy 4J. Accounting for the fi- 
nal state interaction leads to partial renormalization of 
the 2M spectral weight (5|-||] and shifts the maximum of 
the 2M spectral weight down to ^ 2.8J. As a result, 
the observed 2M line shape (See Fig. |l|) consists of two 
peaks: one is close to 4 J, and the other, due to magnon- 
magnon scattering, is at about 2.8 J {Jscoc ~ 130 meV, 
Jybco ~ 110 meV) The relative intensity of these 
two peaks depends on the triple resonance conditions. 

(c) In the Aig channel the long wavelength magnons 
from the vicinity of the magnetic Brillouin zone cen- 
ter contribute to 2M scattering and determine the low- 
frequency slope of the 2M lineshape for Sr2Cu02Cl2. 
The threshold at about 1750 cm^^ in the 2M intensity ob- 
served for YBa2Cu306+5 could be the result of interlayer 
superexchange J_i_ within the Cu02-plane bilayers: the 
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interlayer coupling results in an optical magnon branch 
with a gap of 2\/ JJ± ||3^] and an onset for 2M intensity at 
Ay/Jir[ iQ. The latter gives an estimate for w J/4. 

(d) The triple resonance theory predicts two peaks for 
the RREP at w = 2.8J: the first one at = ljI""'^ « 
2A + 2. 9 J, and the second at ojI'^'^'^ , almost at the upper 
edge of the electron-hole excitation band. This predic- 
tion correlates with observed features (See Fig. where 
the first resonance is about 2.9 — 3.4 J above 2A, and the 
RREP shows a strong rise to the second resonance. More- 
over, the model calculation of the Raman vertex in the 
absence of damping predicts that the 2M peak intensity 
increases by an inverse linear law as the incoming photon 
energy approaches the upper edge of the charge-transfer 
excitation band. In the inset of the Fig. || we display the 
satisfactory linear fit of the inverse 2M intensity data for 
YBa2Cu306+5 versus Wi and extract the linear singular- 
ity excitation energy « 3.1 eV. The latter allows 
one to estimate the electron-hole excitation band width 
as a;['^^^ — 2 A « 1.3 eV and nearest-neighbor hopping 
integral tyBCO ~ 320 meV since ul"'^ « 2VA2 + 16*2. 

Despite the agreement with experimental results, the 
triple resonance calculations |^,^ use a simplified effec- 
tive one-band Hubbard model which takes into account 
only nearest-neighbor hopping integrals t. By breaking 
electron-hole symmetry, the Hubbard model with next- 
nearest-neighbour hopping t' |36| provides satisfactory 
agreement between the narrow 0.3 eV) valence band 
width observed in photoemission experiments [ p7[ and 
our estimate 1.3 eV) for the sum of the valence and 
conduction bands width. Accounting for t' also results in 
a different strength and position of the triple resonance 
at ujI^'^^ in the two investigated AF's, but does not effect 
ujI'^'^^ and our estimate of the hopping integral t . 

In summary, we have performed low-temperature 
two-magnon resonance Raman excitation profile mea- 
surements for single layer Sr2Cu02Cl2 and bilayer 
YBa2Cu306-i-5 AF insulators for the broad excitation re- 
gion from 1.65 to 3.05 eV. We have shown that the ex- 
perimental data are in general agreement with the the- 
ory proposed by Chubukov and Frenkel [p|j27| , indicating 
that the triple resonance process gives the leading con- 
tribution to resonant 2M scattering. This process con- 
tributes to Big, Aig and i?2g scattering geometries. The 
Big 2M line shape is composite and contains two peaks 
at about 2.8 and 4J. The intensities of the two peaks 
strongly depend on the excitation resonance conditions. 
By measuring the excitation energy of the triple reso- 
nance ojI^'^'^ and knowing the CTG from the optical con- 
ductivity, we estimate the electron-hole excitation band 
width to be « 1.3 eV and the nearest-neighbor hopping 
integral tyBCO ~ 320 meV. The Aig 2M band of bilayer 
YBa2Cu306-i-5 material exhibits an optical gap that gives 
an estimate of interlayer superexchange J± w J/4. 
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FIG. 1. The two-magnon resonance Raman scattering 
cross sections {x'y' and x'x' geometries) for Sr2Cu02Cl2 and 
YBa2Cu306+4 single crystals for different excitation ener- 
gies above the charge transfer gap 2A. The scattering in- 
tensities for Sr2Cu02Cl2 are reduced by factor of 10 rela- 
tive to YBa2Cu306+4. Baselines are shifted as indicated. 
Dashed lines mark 2.8J -I- J± and 4J + J± Raman shifts. 
Jscoc = 1060 cm~\ Jybco = 895 cm~\ J±scoc = 0, 
J±YBCO = 215 cm"\ 

FIG. 2. The intensity of the Big two-magnon peak at 2.8 J 
(circles) as a function of the excitation energy uji. Also shown 
is the imaginary part of the dielectric constant from Refs. 
[ po|j3H (solid line). Dashed lines mark the charge transfer gap 
energy 2 A and 2A -I- 3 J. Arrows indicate the triple resonance 
excitation energies. The inset shows the linear fit (dashed 
line) to the inverse two-magnon peak intensity (open circles). 
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